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Abstract. The non-linear dynamics of self-gravitating irrotational dust 
is analyzed in a general relativistic framework, using synchronous and co- 
moving coordinates. Writing the equations in terms of the metric tensor 
of the spatial sections orthogonal to the fluid flow allows an unambiguous 
expansion in inverse powers of the speed of light. The Newtonian and 
post-Newtonian approximations are derived in Lagrangian form. A gen- 
eral formula for the gravitational waves generated by the non-linear evo- 
lution of cosmological perturbations is given. It is argued that a stochastic 
gravitational-wave background is produced by non-linear cosmic struc- 
tures, with present-day closure density Q gw ~ 1CT 5 - 10~ 6 on Mpc scale. 



1. Introduction 

The gravitational instability of collisionless matter in a cosmological framework 
is usually studied within the Newtonian approximation, which basically con- 
sists in neglecting terms higher than the first in metric perturbations around 
a matter-dominated Friedmann-Robertson- Walker (FRW) background, while 
keeping non-linear density and velocity perturbations. This approximation is 
usually thought to produce accurate results in a wide spectrum of cosmological 
scales, namely on scales much larger than the Schwarzschild radius of collapsing 
bodies and much smaller than the Hubble horizon scale, where the peculiar grav- 
itational potential ip g , divided by the square of the speed of light c 2 to obtain a di- 
mensionless quantity, keeps much less than unity, while the peculiar matter flow 
never becomes relativistic. To be more specific, the Newtonian approximation 
consists in perturbing only the time-time component of the FRW metric tensor 
by an amount 2ip g /c 2 , where ip g is related to the matter density fluctuation 5 via 
the cosmological Poisson equation, V 2 ip g (x, r) = 4irGa 2 (T)gb{T)5(x, r), where 
Qb is the background matter density, a(r) the appropriate FRW scale-factor and 
r the conformal time. The fluid dynamics is then usually studied in Eulerian 
coordinates by accounting for mass conservation and using the cosmological ver- 
sion of the Euler equation for a self-gravitating pressureless fluid to close the 
system. To motivate the use of this "hybrid approximation", which deals with 
perturbations of the matter and the geometry at a different perturbative or- 
der, one can either formally expand the correct equations of General Relativity 
(GR) in inverse powers of the speed of light or simply notice that the peculiar 
gravitational potential is strongly suppressed with respect to the matter per- 
turbation by the square of the ratio of the perturbation scale A to the Hubble 
radius r# = off -1 (H being the Hubble constant): (p g /c 2 ~ 5 (X/rn) 2 - 
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Such a simplified approach, however, already fails in producing an accurate 
description of the trajectories of relativistic particles, such as photons. Neglect- 
ing the relativistic perturbation of the space-space components of the metric, 
which in the so-called longitudinal gauge is just —2(p g /c 2 , would imply a mistake 
by a factor of two in well-known effects such as the Sachs-Wolfe, Rees-Sciama 
and gravitational lensing. The level of accuracy not only depends on the pecu- 
liar velocity of the matter producing the spacetime curvature, but also on the 
nature of the particles carrying the signal to the observer. Said this way, it 
may appear that the only relativistic correction required to the usual Eulerian 
Newtonian picture is that of writing the metric tensor in the "weak field" form 
(e.g. Peebles 1993) 



As we are going to show, this is not the whole story. It is well-known 
in fact that the gravitational instability of aspherical perturbations (which is 
the generic case) leads to the formation of very anisotropic structures whenever 
pressure gradients can be neglected (e.g. Shandarin et al. 1995 and references 
therein). Matter first flows in almost two-dimensional structures called pan- 
cakes, which then merge and fragment to eventually form one-dimensional fil- 
aments and point-like clumps. During the process of pancake formation the 
matter density, the shear and the tidal field formally become infinite along 
evanescent two-dimensional configurations corresponding to caustics; after this 
event a number of highly non-linear phenomena, such as vorticity generation by 
multi-streaming, merging, tidal disruption and fragmentation, occur. Most of 
the patology of the caustic formation process, such as the local divergence of the 
density, shear and tide, and the formation of multi-stream regions, are just an 
artifact of extrapolating the pressureless fluid approximation beyond the point 
at which pressure gradients and viscosity become important. In spite of these 
limitations, however, it is generally believed that the general anisotropy of the 
collapse configurations, either pancakes or filaments, is a generic feature of cos- 
mological structures originated through gravitational instability, which would 
survive even in the presence of a collisional component. 

This simple observation shows the inadequacy of the standard Newtonian 
paradigm. According to it the lowest scale at which the approximation can be 
reasonably applied is set by the amplitude of the gravitational potential and 
is given by the Schwarzschild radius of the collapsing body, which is negligibly 
small for any relevant cosmological mass scale. What is completely missing in 
this criterion is the role of the shear which causes the presence of non-scalar 
contributions to the metric perturbations. A non-vanishing shear component 
is in fact an unavoidable feature of realistic cosmological perturbations and af- 
fects the dynamics in (at least) three ways, all related to non-local effects, 
i.e. to the interaction of a given fluid element with the environment. First, 
at the lowest perturbative order the shear is related to the tidal field gener- 
ated by the surrounding material by a simple proportionality law. Second, it 
is related to a dynamical tidal induction: the modification of the environment 
forces the fluid element to modify its shape and density. In Newtonian gravity, 
this is an action-at-a-distance effect, which starts to manifest itself in second- 
order perturbation theory as an inverse-Laplacian contribution to the velocity 
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potential (e.g. Catelan et al. 1995). Third, and most important here, a non- 
vanishing shear field leads to the generation of a traceless and divergenceless 
metric perturbation which can be understood as gravitational radiation emitted 
by non-linear perturbations. This contribution to the metric perturbations is 
statistically small on cosmologically interesting scales, but it becomes relevant 
whenever anisotropic (with the only exception of exactly one-dimensional) col- 
lapse takes place. In the Lagrangian picture such an effect already arises at the 
post-Newtonian (PN) level. Note that the two latter effects are only detected 
if one allows for non-scalar perturbations in physical quantities. Contrary to a 
widespread belief, in fact, the choice of scalar perturbations in the initial condi- 
tions is not enough to prevent tensor modes to arise beyond the linear regime in 
a GR treatment. Truly tensor perturbations are dynamically generated by the 
gravitational instability of initially scalar perturbations, independently of the 
initial presence of gravitational waves. This point is very clearly displayed in 
the GR Lagrangian second-order perturbative approach. The pioneering work 
in this field is by Tomita (1967), who calculated the gravitational waves ir a g 
emitted by non-linearly evolving scalar perturbations in an Einstein-de Sitter 
background, in the synchronous gauge. Matarrese, Pantano & Saez (1994a,b) 
obtained an equivalent result but with a different formalism in comoving and 
synchronous coordinates. 

Recently a number of different approaches to relativistic effects in the non- 
linear dynamics of cosmological perturbations have been proposed. Matarrese, 
Pantano & Saez (1993) proposed an algorithm based on neglecting the magnetic 
part of the Weyl tensor in the dynamics, obtaining strictly local fluid-flow evo- 
lution equations, i.e. the so-called "silent universe". This formalism, however, 
cannot be applied to cosmological structure formation inside the horizon, where 
the non-local tidal induction cannot be neglected, i.e. the magnetic Weyl ten- 
sor H a p is non-zero, with the exception of highly specific initial configurations 
(Matarrese et al. 1994a; Bertschinger & Jain 1994; Bruni, Matarrese & Pantano 
1995a; the dynamical role of H a g was also discussed by Bertschinger & Hamilton 
1994 and Kofman & Pogosyan 1995). Rather, it is probably related to the non- 
linear dynamics of an irrotational fluid outside the (local) horizon (Matarrese 
et al. 1994a, b). One possible application (Bruni, Matarrese & Pantano 1995b), 
is in fact connected to the Cosmic No-hair Theorem. Matarrese & Terranova 
(1995) followed the more "conservative" approach of expanding the Einstein 
and continuity equations in inverse powers of the speed of light, which then 
defines a Newtonian limit and, at the next order, post-Newtonian corrections. 
Their approach differs from previous ones, because of the gauge choice: we used 
synchronous and comoving coordinates, because of which this approach can be 
called a Lagrangian one. Various approaches have been proposed in the liter- 
ature, which are somehow related. A PN approximation has been followed by 
Futamase (1991) to describe the dynamics of a clumpy universe. Tomita (1991) 
used non-comoving coordinates in a PN approach to cosmological perturbations. 
Shibata & Asada (1995) recently developed a PN approach to cosmological per- 
turbations, also using non-comoving coordinates. Kasai (1995) analyzed the 
non-linear dynamics of dust in the synchronous and comoving gauge. 
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2. Method 

We consider a pressureless fluid with vanishing vorticity. Using synchronous and 
comoving coordinates, the line-element reads 

ds 2 = a 2 (r)[- c 2 dr 2 + la/3 (q, r)dq a dq^ , (2) 

where we have factored out the scale-factor of the isotropic FRW solutions. 

By subtracting the isotropic Hubble-flow, we introduce a peculiar velocity- 
gradient tensor -& a p = ^7 a7 7 7/ g', where primes denote differentiation with respect 

tO T. 

Thanks to the introduction of this tensor we can write the Einstein's equa- 
tions in a cosmologically convenient form. The energy constraint reads 

$ 2 - + + c 2 {lZ - 6k) = 167rGa 2 g b 6 , (3) 

where lZ a p(pf) is the conformal Ricci curvature of the three-space with met- 
ric 7 a/ 3; for the background FRW solution J^B^ = (1 + § Q 2 )~ 2 ^a3, one has 
lZ a p{^ FRW ) = 2K5 a p. We also introduced the density contrast 5 = (g — Qb)/ Qb- 
The momentum constraint reads 

= ■ ( 4 ) 

The double vertical bars denote covariant derivatives in the three-space with 
metric j a/3 . 

Finally, after replacing the density from the energy constraint and sub- 
tracting the background contribution, the extrinsic curvature evolution equation 
becomes 



ATl a f} -{n + 2K)8 a p 



0. (5) 



The Raychaudhuri equation for the evolution of the peculiar volume-expansion 
scalar i? reads 



<?' + ^-tf + + AitGa 2 eb 5 = . (6) 

The main advantage of this formalism is that there is only one dimensionless 
(tensor) variable in the equations, namely the spatial metric tensor j a p. The 
only remaining variable is the density contrast which can be written in the form 

%» = (1 + 5 (q))h(q,T)ho(q)]- 1/2 - 1 , (7) 
where 7 = det 7 Q( g. 



3. Results and conclusions 

The method is then based on a 1/c 2 expansion of equations above which first of 
all leads to a new, purely Lagrangian, derivation of the Newtonian approximation 
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(Matarrese &; Terranova 1995). One of the most important result in this respect 
is that we obtained a simple expression for the Lagrangian metric; exploiting the 
vanishing of the spatial curvature in the Newtonian limit we were able to write 
it in terms of the displacement vector S(q, r) = x(q, r) — q, from the Lagrangian 
coordinate q to the Eulerian one x of each fluid element (e.g. Buchert 1995 and 
references therein), namely 



ds 2 = a 2 { T ) 



-cdr +8 AB [8 a + —^)[S fi + —^ r - 



(8) 



A straightforward application of this formula is related to the Zel'dovich approx- 
imation. The spatial metric is that of Euclidean space in time-dependent curvi- 
linear coordinates, consistently with the intuitive notion of Lagrangian picture 
in the Newtonian limit. Read this way, the complicated equations of Newto- 
nian gravity in the Lagrangian picture become much easier: one just has to deal 
with the spatial metric tensor and its derivatives. The displacement vector is 
then completely fixed by solving the Raychaudhuri equation together with the 
momentum constraint in the c — > oo limit. 

Next, we can consider the post-Newtonian corrections to the metric and 
write equations for them. In particular, we can derive a simple and general equa- 
tion for the gravitational-waves 7r Q/ g emitted by non-linear structures described 
through Newtonian gravity. The result can be expressed both in Lagrangian 
and Eulerian coordinates. In the latter case one has, 



an coordinates. In the latter case one has, 
V 2 x ir AB = $1% + 6abV 2 x 9W + 2(Mab ~ $ac$ C b) , 



(9) 



with capital latin labels A,B,... = 1,2,3 indicating Eulerian coordinates and 
V^ii^ = — \{"& 2 — 'd A B 'd B A ), which generally allows a simple derivation of ttab, 
given the (gradients of the) velocity potential, &ab = d 2 <& v /dx A dx B , by a con- 
volution in Fourier space. These formulae would allow to calculate the amplitude 
of the gravitational-wave modes in terms of the velocity potential, which in turn 
can be deduced from observational data on radial peculiar velocities of galaxies. 

In the standard case, where the cosmological perturbations form a homoge- 
neous and isotropic random field, we can obtain a heuristic perturbative estimate 
of their amplitude in terms of the rms density contrast and of the ratio of the 
typical perturbation scale A to the Hubble radius r# = off -1 . One simply 
has TT rrns /c 2 ~ 5 2 ms (X/rH ) 2 . This effect gives rise to a stochastic background 
of gravitational waves which gets a non-negligible amplitude in the so-called 
extremely -low-frequency band (e.g. Thorne 1995), around 10~ 14 - 10~ 15 Hz. We 
can roughly estimate that the present-day closure density of this gravitational- 
wave background is 

n gw (\) ~ C s (^) 2 • (io) 

In standard scenarios for the formation of structure in the universe, the typical 
density contrast on scales 1-10 Mpc implies that Cl gw is about 10 -5 - 1CT 6 . 
We might speculate that such a background would give rise to secondary CMB 
anisotropies on intermediate angular scales: a sort of tensor Rees-Sciama effect. 
This issue will be considered in more detail elsewhere. 
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The previous PN formula also applies to isolated structures, where the den- 
sity contrast can be much higher than the rms value, and shear anisotropies 
play a fundamental role. A calculation of Tr a p in the case of a homogeneous 
ellipsoid showed that the PN tensor modes become dominant, compared to the 
Newtonian contributions to the metric tensor, during the late stages of collapse, 
and possibly even in a shell-crossing singularity. It is important to stress that 
this effect generally contradicts the standard paradigm that the smallest scale 
for the applicability of the Newtonian approximation is set by the Schwarzschild 
radius of the object. Such a critical scale is indeed only relevant for nearly spher- 
ical collapse, whereas this effect becomes important if the collapsing structure 
strongly deviates from sphericity. 
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